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1.  Introduction 


In  support  of  Kiowa  Warrior  (KW)  main  rotor  blade  ballistic  vulnerability  tests  under  the  Live 
Fire  Test  and  Evaluation  program  (7),  a  comprehensive  helicopter  aeroelastic  analysis  was 
conducted  to  calculate  flight  loads  on  the  OH-58D  KW  main  rotor  blade  at  hover  and  forward 
flight  conditions.  Centrifugal  force  (CF)  and  flap  and  chord  bending  moments  were  calculated 
with  the  University  of  Maryland  Advanced  Rotorcraft  Code  (UMARC)  comprehensive 
helicopter  aeroelastic  analysis  code  (2). 

The  calculated  loads  were  first  correlated  with  available  flight  data  for  the  Bell  Model  406LM 
helicopter  to  validate  the  analytical  model.  After  the  correlation  study,  UMARC  was  used  to 
calculate  blade  loads  at  the  ballistic  test  section  radial  stations.  These  baseline  loads  were 
represented  during  the  ballistic  tests  (7). 

The  UMARC  helicopter  aeroelastic  model  used  in  the  present  study  is  a  nonlinear  representation 
of  elastic  rotor  blades  coupled  to  a  rigid  fuselage.  The  key  feature  in  UMARC  is  the  use  of  a 
finite  element  (FE)  methodology,  which  can  accurately  model  the  kinematics  and  elastic 
behavior  of  the  rotating  blade.  The  rotor  blade  is  assumed  to  be  an  elastic  beam  undergoing  flap 
bending,  lag  bending,  elastic  twist,  and  axial  deflections.  The  helicopter  fuselage  is  assumed  to 
be  a  rigid  body  undergoing  six  degrees  of  freedom.  Each  blade  is  modeled  as  a  number  of  beam 
EEs.  The  formulation  for  the  blade  and  fuselage  equations  of  motion  is  based  on  Hamilton’s 
energy  principle.  The  analysis  is  developed  for  helicopters  with  nonuniform  rotor  blades  having 
pretwist,  precone,  and  chord-wise  offsets  of  the  center  of  mass,  aerodynamic  center,  and  tension 
center,  from  the  elastic  axis. 

This  report  presents  the  methods  and  results  of  the  aeroelastic  analysis.  Definitions  of  symbols 
used  in  the  report  are  given  in  the  Fist  of  Abbreviations  and  Symbols.  In  section  2,  the 
formulation  details  and  aerodynamic  model  used  in  UMARC  helicopter  aeroelastic  analysis  code 
are  presented.  Numerical  results  are  shown  and  compared  with  available  flight  data  in  section  3. 
The  summary  and  recommendations  for  future  work  are  presented  in  section  4. 


2.  Helicopter  Aeroelastic  Model 


The  calculation  of  rotor  loads  is  a  difficult  task  because  of  the  complex  interactions  of  structural, 
inertial,  and  aerodynamic  forces  acting  on  the  rotary  wing.  Basically,  it  is  necessary  to  calculate 
the  periodic  aerodynamic  and  inertial  forces  of  the  blade,  and  thus  the  resulting  motion  of 
structural  components.  Since  the  higher  harmonic  blade  loading  is  the  principal  source  of  high 
loads,  an  accurate  analysis  of  the  rotor  aerodynamics  is  required,  including  the  effects  of  the 
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rotor  wake,  stall,  and  compressibility.  The  high  frequencies  involved  and  the  importance  of 
resonance  excitation  also  require  good  inertial  and  structural  models. 

2,1  UMARC  Helicopter  Aeroelastic  Analysis  Code 

The  UMARC  helicopter  aeroelastic  model  used  in  this  study  is  a  nonlinear  representation  of 
elastic  rotor  blades  coupled  to  a  rigid  helicopter  fuselage.  The  key  feature  in  UMARC  is  the  use 
of  an  FE  methodology,  which  can  accurately  model  the  kinematics  and  elastic  behavior  of  the 
rotating  blade.  The  rotor  blade  is  assumed  to  be  an  elastic  beam  undergoing  flap  bending,  lag 
bending,  elastic  twist,  and  axial  deflections.  The  helicopter  fuselage  is  assumed  to  be  a  rigid 
body  undergoing  six  degrees  of  freedom.  Each  blade  is  modeled  as  a  number  of  beam  EEs.  The 
formulation  for  the  blade  and  fuselage  equations  of  motion  is  based  on  Hamilton's  energy 
principle.  The  analysis  is  developed  for  helicopters  with  nonuniform  rotor  blades  having 
pretwist,  precone,  and  chord-wise  offsets  of  the  center  of  mass,  aerodynamic  center,  and  tension 
center  from  the  elastic  axis.  This  code  is  the  current  state-of-the-art  rotorcraft  engineering 
analysis  code  and  is  widely  used  by  U.S.  Army  Research  Eaboratory  researchers  and  other 
Govemment/industry  researchers  (e.g.,  NASA  Ames  and  Boeing  and  Sikorsky). 

In  UMARC,  aerodynamic  loads  acting  on  the  rotor  blade  are  calculated  using  quasi-steady  strip 
aerodynamic  theory.  Noncirculatory  aerodynamic  forces  are  also  included.  To  include  the  effect 
of  high  angle-of-attack  flows  on  the  retreating  blade,  a  dynamic  stall  model  is  incorporated. 
Dynamic  stall  characterizes  the  delay  in  flow  separation  because  of  unsteady  angle  of  attack  and 
the  shedding  of  a  vortex  from  the  leading  edge  of  the  airfoil  when  it  gets  into  a  deep  stall 
condition.  These  effects  are  introduced  in  the  calculation  of  section  lift  and  drag  and  pitching 
moment.  The  first  step  to  calculate  blade  loads  and  moment  is  to  determine  the  trim  state  of  the 
helicopter  for  given  flight  conditions  (i.e.,  gross  weight,  ambient  conditions,  and  forward  speed). 
The  UMARC  utilizes  the  coupled  trim  solution  procedure  to  calculate  the  helicopter  trim  in 
hover  and  forward  flight. 

The  coupled  trim  analysis  consists  of  two  phases,  vehicle  trim  and  steady  response,  calculated  as 
one  coupled  solution  using  a  modified  Newton’s  method.  Eor  a  given  flight  condition,  the 
control  settings  and  the  blade  steady  response  must  satisfy  both  the  blade  and  the  vehicle 
equilibrium  conditions.  The  method  of  solving  simultaneously  the  blade  responses  and  the  trim 
control  settings  is  referred  to  as  the  coupled  trim  analysis.  An  uncoupled  trim  solution  based  on 
the  rigid  flapping  blade  assumption  is  used  as  an  initial  estimate  for  the  coupled  trim  analysis. 
With  the  trim  control  settings,  the  blade  steady  responses  are  calculated.  Using  the  blade 
responses,  the  rotor  hub  loads  and  a  new  vehicle  equilibrium  position  are  recomputed.  The 
control  settings  are  then  updated  based  on  the  new  equilibrium  condition  (5). 
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2.2  Vehicle  Trim 


Propulsive  trim,  which  simulates  an  aircraft  free-flight  condition,  is  used  to  calculate  the  initial 
rotor  control  settings.  The  solution  is  determined  from  the  overall  equilibrium — three  force 
(vertical,  longitudinal,  and  lateral)  and  three  moment  (pitch,  roll,  and  yaw)  equations.  (See  Kim 
[4]  for  more  details.) 


For  a  specified  helicopter  weight  coefficient,  Cw,  and  a  forward  speed,  p,  the  unknown  quantities 
to  be  determined  from  the  vehicle  equilibrium  equations  are  as  follows: 


\f^  s  s  \c  \s  tail  ]  ' 


(1) 


These  values  are  recalculated  iteratively  using  the  modified  hub  forces  and  moments,  including 
the  blade  elastic  responses.  The  solution  technique  is  based  on  a  modified  Newton’s  method. 
The  rotor  controls,  which  are  updated  at  the  rth  iteration,  can  be  expressed  in  a  canonical  form  as 
follows: 


^  u.  +  Am,  ,  (2) 

in  which 

A»,  F(a,),  (3) 

ou 

where  Uo  is  the  trim  control  settings  obtained  initially  using  the  rigid  flapping  blade  solution  and 
dF 

the  Jacobian  —  is  calculated  using  the  finite  difference  approach.  For  computational 
du 

efficiency,  the  Jacobian  is  computed  only  once  initially  and  used  for  subsequent  iterations. 

2,3  Rotor  Steady  Response 

The  rotor  dynamic  response  involves  the  determination  of  time-dependent  blade  positions  at 
different  azimuth  locations  for  one  rotor  revolution.  To  reduce  computational  time,  the  FE 
equations  are  transformed  into  modal  space  as  a  few  normal  mode  equations  using  the  coupled 
natural  vibration  characteristics  of  the  blade.  These  nonlinear  periodic  coupled  equations  are 
solved  for  steady  response  using  an  FE  in  time  procedure  based  on  Hamilton’s  principle  in  weak 
form. 

One  rotor  revolution  is  divided  into  a  number  of  azimuthal  (time)  elements,  and  then  periodicity 
of  response  is  used  to  join  the  motions  of  the  first  and  last  elements.  The  assembly  of  elements 
results  in  nonlinear  algebraic  equations  that  are  solved  using  the  Newton-Raphson  procedure. 

After  the  blade  response  is  obtained,  hub  forces  and  moments  are  calculated  using  the  force 
summation  method.  These  values  are  finally  used  to  recalculate  the  coupled  trim  control  values, 
as  described  in  section  2.2. 
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2.4  Rotor  Wake  Modeling 

For  the  induced  in-flow  distribution  on  the  rotor  disk,  a  free-wake  model  is  fully  coupled  in  the 
rotor  aeroelastic  analysis.  The  model  can  account  for  wake  self-distortion  by  updating  its 
geometry  according  to  newly  calculated  in-flow  and  blade  circulation  distributions.  The 
geometry  of  the  free  wake  is  divided  into  the  following  three  regions:  (1)  near  wake,  (2)  rolling 
up  wake,  and  (3)  far  wake.  The  near  wake  consists  of  a  series  of  radial  panels,  each  with  linear 
circulation  distributions.  The  rolling-up  wake  consists  of  an  inboard  linear  circulation 
distribution  panel,  and  a  tip  panel  that  represents  the  rolling  up  of  the  tip  vortex.  The  far  wake  is 
modeled  as  one  panel  with  a  linear  circulation  distribution  and  a  concentrated  tip  vortex  whose 
strength  is  proportional  to  the  maximum  circulation  value  on  the  rotor  blade.  The  helical 
geometry  of  the  concentrated  tip  vortex  is  updated,  while  the  inboard  wake  portions  are  not 
changed. 

The  free-wake  analysis  is  implemented  in  three  stages.  First,  blade  motion  and  loading  are 
calculated  using  a  linear  in-flow  model.  Next,  wake-induced  coefficients  are  calculated  for 
undistorted  wake  geometry.  The  nonuniform  inflow  is  calculated  and  used  to  obtain  blade 
motion  and  loading.  Finally,  the  free-wake  geometry  is  calculated.  For  this,  the  influence 
coefficients  are  reevaluated,  and  blade  motion  and  loading  are  again  obtained  using  nonuniform 
in-flow  values.  For  subsequent  iterations,  the  free-wake  geometry  is  generally  held  fixed,  and 
only  the  strength  of  vortices  is  updated.  This  comprehensive  rotor  wake  model  is  used  in  the 
present  blade  flight  loads  calculation. 

2.5  Solution  Procedure  for  Blade  Load  Calculations 

The  following  procedures  were  used  to  calculate  blade  loads  and  moments  on  selected  blade 
stations  in  hover  and  forward  flight  conditions. 

1 .  With  prescribed  input  data,  vehicle  trim  equations  are  calculated  using  rigid  blade  flapping 
(a  starting  point). 

2.  Using  control  inputs  from  the  vehicle  trim  solution  of  step  1,  the  blade  nonlinear  steady 
response  is  calculated.  The  results  give  detailed  individual  blade  responses  at  different 
span-wise  and  azimuthal  positions. 

3.  Hub  loads  and  moments  are  calculated  using  elastic  rotor  responses.  Then,  the  vehicle  trim 
values  and  blade  responses  are  recalculated  iteratively  using  the  modified  hub  forces  and 
moments.  This  step  is  repeated  until  a  converged  solution  is  obtained. 

4.  Blade  loads  and  moments  at  selected  blade  stations  are  calculated  using  the  force 
summation  method. 
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3.  Results  and  Discussion 


Numerical  results  are  first  calculated  for  the  OH-58D  KW  helicopter  at  the  gross  weight  of 
4500  lb  with  the  forward  speed  of  120  kn  (equivalent  advance  ratio,  q,  of  0.28).  This  flight 
condition  is  selected  for  the  benchmark  correlation  with  the  available  test  data  obtained  on  the 
Bell  Model  406T  rotor  blade  (5,  6).  The  406T  rotor  blade  is  also  used  on  the  OH-58D  KW. 
Following  the  correlation,  blade  loads  were  calculated  for  the  KW  at  a  gross  weight  of  4086  lb 
for  hover  and  a  cruising  forward  speed  of  100  kn  (q  =  0.23). 

The  KW  helicopter  is  equipped  with  a  four-bladed,  soft  in-plane,  advanced  geometry,  hingeless 
main  rotor  system.  Some  important  aerodynamic  characteristics  of  this  helicopter  rotor  are  given 
in  table  1 .  In  the  present  study,  an  effort  was  made  to  model  the  helicopter  main  rotor  blade 
characteristics  as  accurately  as  possible.  The  rotor  blade  structural  properties  data  obtained  from 
manufacture  (7,  8)  were  used  to  create  a  finite  element  model  (FEM)  (see  table  2).  Figure  1 
shows  the  global  and  local  degrees  of  freedom  of  the  FE  blade  model  for  the  whole  rotor  blade. 
Nonlinear  twist  distribution  of  the  blade  was  also  modeled  in  the  blade  EEs. 


Table  1.  Aerodynamic  characteristics  of  KW  helicopter. 


Air  Density, 

P 

(slug/fF) 

No.  of 
Blades 
in 

Main 

Rotor, 

Nb 

Helleopter 

Gross 

Welght(s), 

W 

(lb) 

Main 

Rotor 

Blade 

Radios, 

R 

(ft) 

Blade 

Tip 

Speed, 

QR 

(ft/s) 

Blade 

Mean 

Chord, 

Cm 

(ft) 

Main 

Rotor 

Solidity 

Ratio, 

a 

Blade 

Twist, 

Otw 

(degree) 

Rotor 

Shaft 

Tilt, 

a, 

(degree)‘‘ 

Blade 

Preeone, 

Pp 

(degree) 

Rotor 

Disk 

Area, 

A 

(ft^) 

.002378(ps,) 

4 

4500,  4086 

17.5 

723.88 

0.786 

0.057 

-12 

5 

2 

962 

“Positive  for  forward  tilt. 


Table  2.  FE  beam  model  input  for  OFt-58D  KW  main  rotor  blade. 


Element  No. 

Flap-Wise 
Bending  Stiffness 

Lag-Wise  Bending 
Stiffness 

/  \ 

El. 

Torsional  Stiffness 

r  Gj  ] 

Element  Length 

UJ 

] 

l(tip) 

0.00581 

0.02316 

0.00522 

0.025 

2 

0.00623 

0.04324 

0.00634 

0.106 

3 

0.00934 

0.06237 

0.00935 

0.169 

4 

0.00823 

0.05762 

0.00843 

0.100 

5 

0.01420 

0.08321 

0.01743 

0.152 

6 

0.02342 

0.16342 

0.02654 

0.162 

7 

0.03429 

0.32143 

0.03021 

0.086 

8  (root) 

0.22840 

0.72342 

0.06443 

0.084 
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For  the  calculation  of  blade  dynamic  response,  each  rotor  blade  is  discretized  into  eight  beam 
elements,  and  each  beam  element  consists  of  15  nodal  degrees  of  freedom  (see  figure  1).  For 
normal  mode  reduction,  five  coupled  rotating  natural  modes  comprised  of  two  flaps,  two  lags, 
and  one  torsion  mode  are  used.  For  periodic  response,  one  cycle  of  time  is  discretized  into  eight 
time  elements,  and  each  time  element  represents  a  fifth-order  polynomial  distribution  of  motion. 


Local  Degrees  of  Freedom  for  15-Node  Finite  Element  Beam  Model 


Global  Degrees  of  Freedom  for  One  Blade  with  8  Beam  Elements 

®  @  ©  ®  @  (D  (2) 


72.71.70 

6.^.62.6I 

54.  55. 52 

4144.4.3 

.36.  .35.  .34 

27. 26. 25 

IS. 17. 16 

y.  K.  7 

64-69  55-60  46-51  37-42  28-33  19-24  10-15 


Figure  1 .  Local  and  global  degrees  of  freedom  of  KW  main  rotor  blade. 

Figure  2  shows  the  correlation  of  the  blade  rotating  natural  frequencies.  The  UMARC  calculated 
values  are  compared  with  the  experimental  data  (7)  obtained  on  the  Bell  406  helicopter  rotor 
blade.  The  UMARC  results  are  calculated  with  the  nominal  operating  rotational  value  of 
395  rpm.  There  is  a  good  correlation  between  the  analytical  UMARC  results  and  the 
experimental  data,  especially  for  the  lower  modes.  It  is  also  seen  that  UMARC  slightly 
underpredicts  the  torsional  frequency. 
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Figure  2.  Correlation  of  main  rotor  blade  natural  frequencies. 

3,1  Correlation  Study 

To  validate  the  current  blade  load  analysis,  the  calculated  blade  centrifugal  loads  and  bending 
moments  were  compared  with  available  flight  test  data.  The  correlation  study  will  not  only  serve 
as  an  analysis  check,  but  also  as  a  guideline  for  further  UMARC  model  improvement  efforts. 

Figure  3  shows  the  correlation  of  the  centrifugal  force  acting  on  the  blade  stations.  The  UMARC 
calculated  values  are  compared  with  the  flight  test  data  obtained  on  the  Bell  406  helicopter.  The 
centrifugal  forces  measured  during  the  flight  test  are  available  at  the  two  blade  stations — 

26,600  lb  at  inboard  blade  station  24.25  (r/R  =  0.1 16)  (5)  and  20,600  lb  at  outboard  blade  station 
96  (r/R  =  0.457)  (6).  Comparing  the  UMARC  results  obtained  using  nonuniform  FE  beam 
model  with  the  flight  test  data,  we  observe  an  excellent  correlation.  Also  shown  in  the  figure  is 
the  analytical  result  obtained  with  a  simple  uniform  blade  assumption.  As  expected,  the  result 
obtained  using  the  uniform  blade  shows  a  poor  correlation.  This  is  most  likely  due  to  the  high 
nonuniformity  of  the  helicopter  blade;  rotor  analysis  with  a  simple  uniform  blade  model  can 
generate  erroneous  results. 
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Correlation  of  Axial  Blade  Loads  on  OH-58D  KW  Blade 


Figure  3.  Correlation  of  CFs  along  the  blade  span. 

The  computed  blade  bending  moments  were  also  correlated  with  the  flight  test  data.  In  reference 
(5),  three  flap  bending  moment  values  are  given  at  three  inboard  blade  stations — stations  3 1 
(r/R  =  0.148),  40  (r/R  =  0.191),  and  50  (r/R  =  0.238),  respectively.  In  reference  (6),  only  one 
bending  value  is  given  at  the  outboard  blade  station  96  (r/R  =  0.457).  These  flight  values  were 
replicated  and  shown  in  table  3.  Also,  shown  in  the  fourth  column  of  the  table  are  the  peak-to- 
peak  values  of  the  moments.  The  peak-to-peak  value  is  a  good  indication  of  loading  on  the 
structural  component.  It  is  seen  that  the  bending  moment  is  largest  inboard  and  decreases 
outboard. 

Figure  4  shows  the  correlation  of  the  UMARC  flap-bending  moments  with  the  flight  test  data. 
There  are  some  discrepancies  between  the  UMARC  results  and  the  flight  test  data.  However, 
considering  the  complexities  associated  with  rotor  loads  and  vibration  analysis,  a  reasonably 
good  agreement  is  observed.  In  particular,  the  general  trend  and  the  peak-to-peak  magnitude  of 
the  moments  are  well  correlated  with  flight  test  data,  especially  at  inboard  blade  stations. 
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Table  3.  Flap  bending  moments  of  Bell  model  406T  helieopter  main  rotor  blade. 


Blade  Station 
(r/R) 

Mean 

Valne 

(in-lbf) 

Oscillatory  Valne 
(in-lbf) 

Peak-to-peak 

Valne 

(in-lbf) 

0.148 

2125 

±2725 

5450 

0.191 

1725 

±2360 

4720 

0.238 

1275 

±2000 

4000 

0.457 

-180 

±780 

1560 

Figure  5  shows  the  eorrelation  of  UMARC  lag  bending  moments  with  flight  test  data.  Again, 
there  are  some  discrepancies  between  the  analytical  results  and  the  test  data.  A  comparison 
reveals  that  the  UMARC  values  are  slightly  greater  than  the  test  data  values.  Flowever,  a 
reasonable  correlation  with  the  test  data  is  observed,  and  the  general  trend  is  especially  well 
correlated. 
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Next,  the  UMARC  caleulated  blade  torsional  moments  were  correlated  with  the  flight  test  data  in 
Figure  6.  There  are  only  two  blade  stations  where  the  flight  torsional  data  are  available — the 
inboard  blade  station  24.25  (r/R  =  0.1 16)  (5)  and  the  outboard  blade  station  96  (r/R  =  0.457)  (6). 
The  inboard  station  (r/R=0.1 16)  is  where  the  main  rotor  blade  is  connected  with  the  pitch  link. 
Again,  there  are  some  discrepancies  between  the  two  results.  Flowever,  a  reasonable  correlation 
is  observed. 

3,2  Blade  Flight  Load  Calculations  for  Hover  and  100-kn  Forward  Speed 

From  the  correlation  study,  it  was  determined  that  the  UMARC  code  adequately  calculated  flight 
loads.  Excellent  agreement  with  test  dada  was  observed  for  the  most  dominant  loading  (i.e.,  the 
axial  load),  giving  a  certain  confidence  in  the  helicopter  rotor  subsystem  model.  Next,  a  new  set 
of  bending  moments  was  calculated  for  the  ballistic  test  flight  conditions — the  gross  weight  of 
4086  lb  at  hover  and  a  cruising  forward  speed  of  100  kn  (p,  =  0.23)  for  three  blade  radial 
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Figure  6.  Correlation  of  torsional  moments  along  blade  span  (p  =  0.28). 

stations.  Since  the  new  flight  conditions  had  a  negligible  influence  on  the  centrifugal  force  on 
the  blade,  the  previous  calculated  values  were  used  for  the  axial  loading.  These  loads  and 
moments  (peak  values)  were  used  as  baseline  loads  during  the  ballistic  test  events.  The  results 
are  summarized  in  table  4. 


Table  4.  Reference  blade  loads  and  moments  for  test  phase  I-B. 


Blade 

Station 

(r/R) 

CF 

(lb) 

FB 

Hover 

GW  =  4086  lb 
(lb-in) 

FB 

lOOkn 

GW  =  4086  lb 
(lb-in) 

FB 

120  kn 

GW  =4500  lb 
(lb-in) 

LB 

Hover 

GW  =  4086  lb 
(lb-in) 

LB 

lOOkn 

GW  =  4086  lb 
(lb-in) 

LB 

120  kn 

GW  =  4500  lb 
(lb-in) 

0.286 

22374 

2020 

2219 

2691 

6860 

7605 

8284 

0.614 

16220 

564 

795 

872 

1416 

1743 

2169 

0.843 

8264 

202 

242 

379 

684 

690 

1112 

Notes:  FB  =  flap  bending  moment,  positive  flapping  up;  LB  =  lag  (chord-wise)  bending,  positive  lagging  back. 
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4.  Summary 


An  aeroelastic  analysis  was  conducted  to  calculate  flight  loads  on  the  OH-58D  KW  main  rotor 
blade  for  hover  and  eruise  forward-flight  eonditions.  Centrifugal  force  and  flap  and  lag  bending 
moments  were  calculated  using  the  UMARC  eode.  The  calculated  loads  were  correlated  with 
flight  test  data  for  the  Bell  model  406LM  helicopter  to  validate  the  analysis.  The  correlation 
showed  good  general  agreement  between  the  test  data  and  analytieal  results.  After  the 
eorrelation  study,  blade  loads  were  calculated  for  the  desired  ballistic  test  conditions. 

In  the  present  study,  the  UMARC  comprehensive  helicopter  aeroelastie  analysis,  based  on  FE 
methodology  in  both  space  and  time  coordinates,  was  condueted  to  analytically  calculate  the 
flight  loads.  In  spite  of  it  being  eomputer-intensive,  its  benefit  must  be  appreciated  as  a  powerful 
analytieal  tool  to  validate  simple  design-oriented  analysis. 
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Advanee  ratio, - 
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Rotor  solidity  ratio,  — ^ — — 

;zR 

Elastie  blade  twist  about  the  elastie  axis 

Eateral  shaft  tilt 

Air  density 

Rotor  rotational  speed 
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